Turbulent-laminar patterns in plane Couette flow 



Dwight Barkley 

Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdon^ 

Laurette S. Tuckcrman 
LIMSI-CNRS, BP 133, 91403 Orsay, FrancR 

Regular patterns of turbulent and laminar fluid motion arise in plane Couette flow near the lowest 
Reynolds number for which turbulence can be sustained. We study these patterns using an extension 
of the minimal flow unit approach to simulations of channel flows pioneered by Jimenez and Moin. 
In our case computational domains are of minimal size in only two directions. The third direction 
is taken to be large. Furthermore, the long direction can be tilted at any prescribed angle to the 
streamwise direction. We report on different patterned states observed as a function of Reynolds 
number, imposed tilt, and length of the long direction. We compare our findings to observations in 
large aspect-ratio experiments. 
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I. INTRODUCTION 

In this chapter we consider plane Couette flow - the flow between two infinite parallel plates moving in opposite 
directions. This flow is characterized by a single non-dimensional parameter, the Reynolds number, defined as Re = 
hll/v, where 2h is the gap between the plates, U is the speed of the plates and v is the kinematic viscosity of the fluid. 
Sec figure ^ For all values of Re, laminar Couette flow uc = yx is a solution of the incompressible Navier-Stokcs 
equations satisfying no-slip boundary conditions at the moving plates. This solution is linearly stable at all values 
of Re. Nevertheless it is not unique. In particular, for Re greater than approximately 325 Q, turbulent states are 
found in experiments and numerical simulations. Our interest is in the flow states found as one decreases Re from 
developed turbulent flows to the lowest limit for which turbulence exists. 

Our work is motivated by the experimental studies of Prigent and coworkers l)200ll2002l2003l2005|) on flow in 
a very large aspect-ratio plane Couette apparatus. Near the minimum Re for which turbulence is sustained, they 
find remarkable, essentially steady, spatially-periodic patterns of turbulent and laminar flow. These patterns emerge 



z 



U 



V 



FIG. 1: Plane Couette geometry. Plates separated by a gap 2h move at speeds ±C/. The coordinate system we use is as shown, 
with y = corresponding to midgap. 
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FIG. 2: Turbulent-laminar pattern at Reynolds number 350. Isosurfaces of streamwise vorticity (u — ±0.5) are shown at one 
instant in time. For clarity the bottom plate is shown in black while the top plate is transparent. The streamwise and spanwise 
extent of the region shown are 60 times the plate separation 2h. 



spontaneously from featureless turbulence as the Reynolds number is decreased. Figure shows such a pattern 
from numerical computations presented in this chapter. Two very striking features of these patterns are their large 
wavelength, compared with the gap between the plates, and the fact that the patterns form at an angle to the 
streamwise direction. 

Fluid flows exhibiting coexisting turbulent and laminar regions have a sig nificant history in fluid dynamics. In the 
mid 1960s a state known as spiral turbulence was first discovered 0, 0, l26j| in counter-rotating Taylor-Couette flow. 
This state consists of a tu rbulent and a laminar region, each with a spiral shape. The experiments of Prigent et 

in a very large aspect-ratio Taylor-Couette system showed that in fact the turbulent and 
laminar regi ons form a period i c patt ern, of which the original observations of Coles and van Atta comprised only one 
wavelength. ICros and Le GaJ l)2002[) discovered large-scale turbulent spirals as well, in the shear flow between a sta- 
tionary and a rotating disk. When converted to comparable quantities, the Reynolds-number thresholds, wavelengths, 
and angles are very similar for all of these turbulent patterned flows. 



II. METHODS 



Our co mputational technique |2| extends the minimal flow unit methodology pioneered by Jimenez and Moin (119911) 
and bv lHamilton et all l)l995j) and so we begin by recalling this approach. Turbulence near transition in plane Couette 
and other channel flows is characterized by the cyclical generation and breakdown of streaks by streamwise-oriented 
vortices. The natural streak spacing in the spanwise direction is about A-bh. In the minimal flow unit approach, 
the smallest laterally periodi c dom ain is sought that can sustain this basic turbulent cycle. For plane Couette flow 
at Re = 400. lHamilton et alJ l|l995() determined this to be approximately (L x , L y , L z ) = (Ah, 2h, 6h). This domain is 
called the minimal flow unit (MFU). The fundamental role of the streaks and streamwise vortices is manifested by the 
fact that the spanwise length of the MFU is near the natural spanwise streak spacing. Figure E{a) shows the MFU in 
streamwisc-spanwise coordinates. 

We extend the MFU computations in two ways. First we tilt the simulation domain in the lateral plane at angle 
9 to the streamwise direction [Figure EJb)]. We use x' and z' for the tilted coordinates. We impose periodic lateral 
boundary conditions on the tilted domain. To respect the spanwise streak spacing while imposing periodic boundary 
conditions in x', the domain satisfies L m i sin 9 ~ Ah for 9 > 0. (For 9 = 0, we require L x > > 6ft..) Secondly, we greatly 
extend one of the dimensions, L z >, past the MFU requirement [Figure[3fc)], in practice between 30h and 220/i, usually 
120/i. 

This approach presents two important advantages, one numerical and the other physical. First, it greatly reduces 
the computational expense of simulating large length-scale turbulent-laminar flows. Our tilted domains need only be 
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FIG. 3: Simulation domains. The wall-normal direction y is not seen; L y = 2ft. The bars represent streamwise vortex pairs 
with a spanwise spacing of 4ft. (The vortices are schematic; these are dynamic features of the actual flow.) (a) MFU domain 
of size 6ft x 4ft. (b) Central portion of a domain [on the same scale as (a)] tilted to the streamwise direction, a, a' and /?, 
P' are pairs of points identified under periodic boundary conditions in x'. (c) Full tilted domain with L x i = 10ft, L z i — 120ft, 
9 = 24°. On this scale the MFU domain, shown for comparison, is small. 




History points 



FIG. 4: Simulation domain. The (x',y) grid is the actual spectral-element mesh used for the case L x / = 10. Only part of the 
z' direction is shown. In practice we use 32 history points in the z' direction. 



long perpendicular to the turbulent bands. In the direction in which the pattern is homogeneous, the domains are of 
minimal size, just large enough to capture the streamwise vortices typical of shear turbulence. Second, the approach 
allows us to impose or restrict the pattern orientation and wavelength. We can thereby investigate these features and 
establish minimal conditions necessary to produce these large-scale patterns. 

We now present some further details of our simulations. We consider the incompressible Navier-Stokes equations 
written in the primed coordinate systems. After nondimensionalizing by the plate speed U and the half gap ft, these 
equations become 

du' 1 

^r + (u'-V')u' = -VV + -^V'V infi, (la) 
ot tie 

V'-u' = in Q, (lb) 

where u'(x', t) is the velocity field and p'(x', t) is the static pressure in the primed coordinate system, and V is used 
to indicate that derivatives are taken with respect to primed coordinates. Q is the computational domain. In these 
coordinates, the no-slip and periodic boundary conditions are 

u'(x',y = ±l,z') = ±(cos0,O,sin0) (2a) 
u'(x' + L x ,,y,z') = u'(x',y,z') (2b) 
u'(x',y,z' + L z ,) = u'(x',y,z') (2c) 

The equations arc simulated using the s pec tral-element (x'-y) - Fourier (z ! ) code Prism [l^. We use a spatial 
resolution consistent with previous studies |lft l27l l28j. Specifically, for a domain with dimensions L x i and L y = 2, 
we use a computational grid with close to L x i elements in the x' direction and 5 elements in the y direction. Within 
each element, we usually use 6th order polynomial expansions for the primitive variables. Figure 01 shows a spectral 
element mesh used for the case of L x > = 10. In the z' direction, a Fourier representation is used and the code is 
parallelized over the Fourier modes. Our typical domain has L z i = 120, which we discretize with 1024 Fourier modes 
or gridpoints. Thus the total spatial resolution we use for the L x i x L y x L z i = 10 x 2 x 120 domain can be expressed 
as N x , x N y x N z , = 61 x 31 x 1024. 
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FIG. 5: Space-time diagram. Kinetic energy E(x' = 0, y = 0, z') at 32 equally spaced points in z' in a domain with L x i x L y x 
L z i = 10 x 2 x 120 with tilt 8 = 24°. The Reynolds number is decremented in discrete steps (right). Three long- lasting and 
well-separated laminar regions emerge spontaneously from uniform turbulence as Re is decreased. 

We shall always use (x, y, z) for the original streamwise, cross-channel, spanwisc coordinates (Figure^. We obtain 
usual streamwise, and spanwise components of velocity and vorticity using u — vl cos 8 + w' sin 8 and w = u' sin 8 — 
w' cos 8, and similarly for vorticity. The kinetic energy reported is the difference between the velocity u and simple 
Couette flow u c , i.e. E = - u c ) 2 + v 2 + w 2 ). 

We have verified the accuracy of our simulations in small domains by comparing to prior simulations [ljj. ^ n 
large domains we have examined mean velocities, Reynolds stresses, and correlations in a turbulent-laminar flow 
at Re = 350 and find that these reproduce experimental results from Taylor-Couette and plane Couette 
|ll| flow. While neither experimental study corresponds exactly to our case, the agreement supports our claim 
that our simulations correctly capture turbulent-laminar states. 

The procedure we use to initiate turbulence is inspired by previous investigations of plane Couette flow in a 
perturbed geometry. We recall that laminar plane Couette flow is linearly stable at all Reynolds numbers. It has been 
found, experimentally Q and numerically 0,13, that the presence of a wire 0] or a ribbon 0, [2f| oriented along 
the spanwise direction causes the flow in the resulting geometry to become linearly unstable to either a steady or a 
turbulent state containing streamwise vortices. We simulate such a flow with a ribbon which is infinitesimal in the 
x' direction, occupies 30% of the cross-channel direction y and spans the entire z' direction. At Re = 500, the effect 
of such a ribbon is to produce a turbulent flow quickly without the need to try different initial conditions. Once the 
turbulent flow produced by the ribbon is simulated for a few hundred time units, the ribbon can be removed and the 
turbulence remains. This procedure is used to initialize turbulent states for the simulations to be described below. 



We begin with one of our first simulations, in a domain tilted at angle 8 = 24°. This angle has been chosen to 
be close to that observed experimentally near pattern onset. The simulation shows the spontaneous formation of 
a turbulent-laminar pattern as the Reynolds number is decreased. We initiated a turbulent flow at Re — 500 by 
perturbing laminar Couette flow with a ribbon as described in Section [H] Time zero in Figure |S] corresponds to the 
removal of the ribbon. The flow is simulated for 500 time units at Re = 500 and the kinetic energy E is measured at 
32 points equally spaced in z' along the line x' = y = in the mid-channel shown in Figure 0] The corresponding 32 
time series are plotted at the corresponding values of z' . At Re = 500, there is no persistent large-scale variation in 
the flow, a state which we describe as uniform turbulence. (This is not the homogeneous or fully developed turbulence 



III. APPEARANCE OF TURBULENT-LAMINAR BANDS 



A. Basic phenomenon 
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that exists at higher Reynolds numbers or in domains without boundaries.) At the end of 500 time units, Re is 
abruptly changed to Re = 450 and the simulation continued for another 500 time units. Then Re is abruptly lowered 
to Re — 425 and the simulation is continued for 1000 time units, etc. as labeled on the right in Figure El 

At Re = 350 we clearly see the spontaneous formation of a pattern. Out of uniform turbulence emerge three regions 
of relatively laminar flow between three regions of turbulent flow. (We will later discuss the degree to which the flow 
is laminar.) While the individual time traces are irregular, the pattern is itself steady and has a clear wavelength of 
40 in the z' direction. This Reynolds number and wavelength arc very close to what is seen in the experiments. 

B. Visualizations 

FigurcElshows visualizations of the flow at the final time in FigurcEl Shown are the kinetic energy, streamwise velocity, 
and streamwise vorticity in the midplane between the plates. The computational domain is repeated periodically to 
tile an extended region of the midplane. The angle of the pattern is dictated by the imposed tilt of the computational 
domain. The wavelength of the pattern is not imposed by the computations other than that it must be commensurate 
with i z /120. The vorticity isosurfaccs of this flow field were shown in Figure [21 Spanwisc and cross-channel velocity 
components show similar banded patterns. 

Clearly visible in the center figure are streamwise streaks typical of shear flows. These streaks have a spanwise 
spacing on the scale of the plate separation but have quite long streamwise extent. We stress how these long streaks 
are realized in our computations. A streak seen in Figure typically passes through several repetitions of the 
computational domain, as a consequence of the imposed periodic boundary conditions. In the single tilted rectangular 
computational domain, a single long streak is actually computed as several adjacent streaks connected via periodic 
boundary conditions. 

Figure0shows the streamwise vorticity and velocity fields between the plates. The two leftmost images correspond 
to the same field as in Figure The streamwise vorticity is well localized in the turbulent regions. Mushrooms of 
high- and low-speed fluid, corresponding to streamwise streaks, can be seen in the turbulent regions of flow. Dark 
velocity contours, corresponding to fluid velocity approximately equal to that of the lower and upper moving plates, 
are seen to reach into the center of the channel in the turbulent regions. In the center of the laminar regions, where 
the flow is relatively quiescent (Figure 0, there is very little streamwise vorticity and the streamwise velocity profile 
is not far from that of laminar Couette flow. In particular, no high- or low-speed fluid reaches into the center of the 
channel in these laminar regions. 

In Figure|Slwc show the mean and rms of the streamwise and spanwise velocity components obtained from averages 
over T = 2000 time units. These results show that the mean flow is maximal at the boundaries separating the 
turbulent and laminar regions while the fluctuations are maximal in the midd le of the turbulent bands. This agrees 
with the experimental observations of of Prigent et al. 1)20011120021 l2003l2005|) Note furt her that the regions of high 
fluctuation have approximately the same rhombic shape as the turbulent regions shown bv lColes and van Attal |)1966|) 
in experiments on Taylor-Couettc flow. Finally, Figure |!|1 shows a space-time plot of streamwise velocity along the 
spanwise line x = y = 0. Specifically, data is taken from reconstructed flows as in Figure Time zero in Figure 
corresponds to the field in Figure Time is ta ken do wnward in this figure to allow for comparison with a similar 
figure from the experimental study by Hegseth l)l996t figure 6) showing the propagation of streaks away from the 
center of turbulent regions. Our results agree quantitatively with those of Hegseth. We find propagation of streaks 
away from the center of the turbulent regions with an average spanwise propagation speed of approximately 0.054 
in units of the plate speed U. Translating from the diffusive time units used by Hegseth, we estimate the average 
spanwise propagation speed of streaks in his data to be approximately 0.060 at Reynolds number 420. This space-time 
plot again shows the extent to which there is some small activity in the regions we refer to as laminar. 

C. Average spectral coefficients 

We have determined a good quantitative diagnostic of the spatial periodicity of a turbulent-laminar pattern. We use 
the same data as that presented in Figure[Sl i.e. velocities at 32 points along the line x' = y = in the midplane along 
the long direction, at each interval of 100 time steps: lOOAf = 1. We take a Fourier transform in z' of the spanwise 
velocity w, yielding w m . We take the modulus \w m \ to eliminate the spatial phase. Finally, we average over a time T 
to obtain (|u) m |). Figure ITU1 shows the evolution of (|iD m |) for wavenumbers m = 3, m = 2, m = 1, and m — during 
one of our simulations (shown below in Figure fTTI which is a continuation of that shown in Figure [SJ). As before, 
the vertical axis corresponds to time, and also to Reynolds number, which was decreased in steps of ARe = 10. We 
average successively over T = 10, T = 30, T = 100, T = 300, and T = 1000 and observe the short-term fluctuations 



FIG. 6: Turbulent-laminar pattern at Re — 
350. The kinetic energy, streamwise velocity, and 
streamwise vorticity are visualized in the y — 
plane, midway between and parallel to the mov- 
ing plates. The computational domain (outlined 
in white, tilted at angle 6 — 24°) is repeated peri- 
odically to tile an extended region in x-z coordi- 
nates. Streamwise streaks, with spanwise separa- 
tion approximately Ah, are visible at the edges of 
the turbulent regions. 



FIG. 7: Turbulent-laminar 
pattern at Re = 350 
viewed between the mov- 
ing plates (x' — plane). 
Left plot shows stream- 
wise vorticity. The other 
three plots show contours of 
streamwise velocity at three 
times separated by 100 time 
units (time increasing left to 
right). The vorticity plot 
and the first velocity plot 
correspond to the field seen 
in Figure |S] 
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FIG. 8: Mean and rms velocity fields for the turbulent-laminar pattern. From left to right: mean streamwise velocity, mean 
spanwise velocity, rms streamwise velocity and rms spanwise velocity. The rms velocities are maximal in the lightest regions. 
Only the central half (30 < z' < 90) of the computational domain is shown. 
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FIG. 9: Space-time plot showing dynamics of the turbulent-laminar pattern. Streamwise velocity is sampled along a spanwise 
cut through the flow field (the line x = y — in the reconstructed flows in Figure |SJ- Time zero corresponds to the time of 
Figure The streaks propagate away from the center of the turbulent regions toward the laminar regions. 

gradually disappear, leaving the long-term features which will be discussed in the next section. We have chosen 
T = 500 as the best compromise between smoothing and preserving the detailed evolution. 

IV. DEPENDENCE ON REYNOLDS NUMBER 

We have investigated in detail the Reynolds-number dependence of the 9 = 24° case. To this end, we have carried out 
two simulations, shown in Figure lTTI In each the Reynolds number is lowered at discrete intervals in time, but following 
a different sequence in the two cases. For each case, we present a space-time diagram of E{x' = 0, y = 0,z',t) at 32 
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FIG. 10: Evolution of (|ty m |), which is an average over time T of the modulus of the Fourier transform in the z' direction of 
32 spanwise velocity samples taken along the line (x' = 0,y = 0). The components with wavenumber m = 3 (solid curve), 
m — 2 (long-dashed curve) , m = 1 (short-dashed curve) and m = (dotted curve) can be used as a quantitative diagnostic of a 
turbulent-laminar pattern. For example, the dominance of the m = 3 component indicates a pattern containing three turbulent 
bands. From left to right, the average is taken over T = 10, T = 30, T = 100, T = 300, and T = 1000. 



values of z' . The Reynolds-number sequence is shown on the right of each diagram and the time (up to T = 59, 000) on 
the left. Each space-time diagram is accompanied by a plot showing the evolution of its average spectral coefficients, 
as defined above. 

Careful observation of Figure [S] already shows a laminar patch beginning to emerge at Re = 390, consistent with 
experimental observations: Prigcnt et al. l)200ll I2002l2003l2005[) observed a turbulent-laminar banded pattern with 
wavelength 46 and angle 25° when they decreased Re below Re = 394. The space-time diagram on the left of 
Figure ITU shows a continuation of this simulation. (Here, the Reynolds numbers intermediate between 500 and 350 
are not shown to reduce crowding.) We see a sequence of different states: uniform turbulence and the three-banded 
turbulent-laminar pattern already seen are succeeded by a two-banded pattern (at Re = 310), then a state containing 
a single localized turbulent band (at Re = 300), and finally laminar Couette flow. These features are reflected in the 
average spectral coefficients. The flow evolves from uniform turbulence (all components of about the same amplitude) 
to intermittent turbulence, to a pattern containing three turbulent bands (dominant m = 3 component) and then two 
turbulent bands (dominant m = 2 component), then a single band (dominant m = 1 and m = 2 components), and 
finally becomes laminar (all components disappear). 

In the simulation on the right, the Reynolds number is decreased more slowly. A state with three bands appears 
at Re = 390. (Although a laminar patch already appears at Re = 400, it is regained by turbulence when Re is 
maintained longer at 400; this is not shown in the figure.) Based on the previous simulation shown on the left, we 
had expected the three turbulent bands to persist through Re = 320. However here, instead, we see a rapid loss of 
two bands, leaving only a single turbulent band. This band moves to the left with a well-defined velocity, emitting 
turbulent spurs toward the right periodically in time. Finally, after a time of T = 36000, one of these spurs succeeds 
in becoming a second turbulent band and the two bands persist without much net motion. It would seem that the loss 
of the second band was premature, and that at Re = 320 one band is insufficient. We then resumed the simulation 
on the left, maintaining Re = 320 for a longer time, and found that two bands resulted in this case as well. Both 
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FIG. 11: Two time series at 8 — 24°. The Reynolds number is lowered in steps, but at different instants to generate the two 
evolutions shown. For each case, we show w at 32 points along the z direction on the right and the spectral components (|ix) m |) 
on the left. Left: uniform turbulence is succeeded by the formation of three bands, then two, then a single band (a localized 
state) and finally by laminar Couette flow. Right: two bands disappear almost simultaneously at Re — 320. The remaining 
band moves toward the left, periodically emitting turbulent spurs, of which one finally becomes a second turbulent band. 



simulations show two bands at Re = 310, one band at Re = 300, and laminar Couette flow at Re = 290. 



A. Three states 



The turbulent-laminar patterned states shown in Figure 1111 arc of three qualitatively different types Q ■ We demon- 
strate this by carrying out three long simulations, at three different Reynolds numbers, that are shown in Figure IT2"1 
In this figure, the energy along the line x' = y = for the 32 points in z 1 has been averaged over windows of length 
T = 500 to yield a value shown by the shading of each space-time rectangle. 

The simulations at Re = 350 and Re = 300 arc carried out by increasing the long direction of our domain, L z i, 
in regular discrete increments of 5 from L z i = 50 to L z i = 140. At Re = 350, a single turbulent band is seen when 
L z i = 50. This band divides into two when L 2 < = 65 and a third band appears when L z < = 130: the periodic pattern 
adjusts to keep the wavelength in the range 35-65. This is close to the wavelength range observed experimentally by 
Prigent et al., which is 46-60. When the same protocol is followed at Re = 300, no additional turbulent bands appear 
as L x t is increased. We call the state at Re = 300 localized and note that turbulent spots are reported near these 
values of Re in the experiments. The small L x i of our computational domain does not permit localization in the x' 
direction; instead localized states must necessarily take the form of bands when visualized in the x-z plane. 

The instantaneous integrated kinetic energy profile 

E = J dx' dy E(x',y,z') 

is plotted at the final time for both cases. For Re = 350, E does not reach zero and the flow does not revert to the 
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FIG. 12: Simulations at Re = 350, Re = 300, and Re = 410 illustrating three qualitatively different regimes. For the simulations 
at Re = 350 and Re = 300, L z i is increased from 50 to 140. The state at Re = 350 is periodic: the turbulent band divides as 
L z i is increased to retain a wavelength near 40. The final kinetic energy profile E(z') is bounded away from zero. The state at 
Re = 300 is localized: a single turbulent band persists, regardless of domain size and E(z') decays exponentially to zero away 
from the band. The simulation at Re = 410 is carried out at L z i = 40. The state is intermittent: laminar regions appear and 
disappear and the average spectral coefficients corresponding to m = 1 (solid) and m = (dashed) oscillate erratically. 

simple Couette solution between the turbulent bands, as could also be seen in the earlier visualizations (Figures 
|5J). In contrast, for Re = 300, E decays to zero exponentially, showing that the flow approaches the simple Couette 
solution away from the turbulent band. In this case, there is truly coexistence between laminar and turbulent flow 
regions. 

The simulation at Re = 410 illustrates another type of behavior. In a domain of length L z > = 40, laminar or, rather, 
weakly-fluctuating regions appear and disappear. The spectral coefficients corresponding to m = 1 (wavelength 40) 
and m = osc illate erratically. Sim ilar states at similar Reynolds numbers are reported experimentally by Prigcnt 
and coworkers l)200lll2002l2003l2005i) , where they are interpreted as resulting from noise-driven competition between 
banded patterns at equal and opposite angles, a feature necessarily absent from our simulations. 



11 





1 


I I 

T 


| 










I 








400 




1 k 




I 

L 


LTT I- 


tie 




1 □ 

3 □ 

□ 


2 


L L 


ic - 

LI 
C 


obO 




GO CD 

3 

3 2 
2 2 


o 

c 

□ 2 
L 

I, 


T 

L 

C 

b 




300 


-T 
- 1 


L L 

C C 

I I 


C 

1 











20 


40 


60 


80 



e 

FIG. 13: Patterns as a function of Reynolds number Re and 9, the angle between the L x i direction of our rectangular 
computational domain and the streamwise direction. The domain is of size L x i x L y x L z i = (4/ sm6) x 2 x 120. For each 
angle, upper and lower limits in Re are shown for each regime. T: uniform turbulence (lower limit in Re). I: intermittent 
turbulence. 3: pattern containing three turbulent-laminar bands, each of approximate wavelength 40. 2: pattern containing 
two bands of approximate wavelength 60. L: pattern containing one turbulent region, possibly localized. C: laminar Couette 
flow (no patterns observed below this Re). Open symbols show experimental observations of Prigent et al. Triangles: patterns 
with wavelength between 46 and 50. Squares: patterns with wavelength between 50 and 60. 

V. DEPENDENCE ON ANGLE 
A. Angle survey 

We have explored the angles with respect to the streamwise direction at which a turbulcnt4aminar pattern may exist. 
The results are plotted in Figure We keep L z , = 120 and L x / = 4/ sin (9. The transition from uniform turbulence 
to laminar Couette flow occurs via intermediate states which occupy a decreasing range of Re as 8 is increased. The 
sequence of states seen for increasing 8 at Re = 350 is qualitatively the same as that for decreasing Re at 8 = 24°: 
uniform turbulence at 6 = 0°, a turbulent-laminar pattern with three bands at 8 = 15° to 8 = 24°, two bands for 
8 = 30° and 8 = 45°, a localized state for 8 = 66°, and laminar Couette flow for 8° > 72. Thus far we have obtained 
patterns for angles between 15° and 66° and the nu mber of bands d ecreases with angle. 

Experimental data from Prig cnt and coworkers l)200ll2003l2005|) is also shown in Figure The wavelengths, 
angles, and Reynolds numbers reported ranged from A 2 < = 46.3 and 8 — 25.3° at Re — 394 to A 2 / = 60.5 and 
8 = 37° at Re = 340. In these ranges of angle and Reynolds number, we observe a similar trend, since our wavelength 
(constrained here to be a divisor of L z / = 120) increases from 40 to 60 as the number of bands decreases from 3 
to 2. Between Re = 325 and Re = 280, experiments showed spots, which may correspond to some of the states we 
have labeled as localized in Figure El At present, we do not systematically distinguished localized states from others 
containing one turbulent region but which may not behave like Figure 1121 The threshold for intermittency is also 
difficult to define and to determine. 

The most striking difference between our computations and the experimental data is that the range of angles over 
which we find periodic turbulent-laminar patterns (from 8 — 15° to at least 8 = 45°) is far greater than that seen in 
the experiment. Patterns with angles outside of the experimental range arc likely to be unstable in a large domain in 
which the angle is unconstrained. 

Our computational technique requires that the size of the domain be increased as 8 decreases according to L x i = 
4/ smd in order to respect the spanwise vortex or streak spacing; see Figure|21 Hence the computational cost increases 
with decreasing 8 and for this reason we have not as yet investigated 8 between 15° and 0°. For 8 exactly 0°, this 
trigonometric constraint is lifted, since the streamwise vortices and streaks would not extend diagonally across the 
rectangular domain, but parallel to its boundaries. As 8 increases, the domain size L x i = 4/ sin 6* decreases, as does 
the computation cost. For 8 between 45° and 90°, for which L X ' is between 5.7 and 4, we reduce the number of 
spectral elements in the x' direction from 10 to 4 (sec Figure ^J. 
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FIG. 14: Evolution for 9 = 90°. Two simulations at Re = 400 are shown, in domains with L x i = L z = 4. The domain shown 
on the right has L z i = L x = 220; that on the left has L z i = L x — 6, close to the minimal flow unit (the two are not shown 
to scale). The large domain supports a transient pattern as an intermediate state between uniform turbulence and laminar 
Couette flow, whereas the turbulence persists in the small domain. 

B. Long streamwise direction 

For 9 = 90°, the domain has a long streamwise direction L z i = L x and a short spanwise direction L x i = L z . FigurclTSl 
shows that, for 9 = 90° and L x = 120, we obtain direct decay from uniform turbulence to laminar Couette flow at 
Re = 385. We have varied L x and show the results in Figure When L x = 220, the turbulence is extinguished at 
Re = 400; a transient pattern of wavelength 110 can be seen. But when L x = 6, we fin d that the turbulence persists 
down to a value of Re ~ 370. We recall that the minimal flow unit was proposed bv lHamilton et all l|l995|) as the 
smallest which can support the streak and streamwise- vortex cycle and maintain turbulence; the flow becomes laminar 
when either of the dimensions are reduced below their MFU values. However, Figure ITT1 shows that turbulence can 
also sometimes be extinguished by increasing L x . Simulations in domains with a long streamwise and a short spanwise 
dimension have also been carried out by Jimenez et al. [l4| with the goal of understanding the role of the streamwise 
dimension, e.g. streak length. 

C. Long spanwise direction 

In the case 9 — 0°, the domain has a long spanwise direction (L z = L z i = 120) and a short streamwise direction 
(L x = L x i = 10). In decreasing the Reynolds number by ARe = 5 after each interval of T = 1000, we observe 
turbulent regions far below Re = 300, terminating only at Re = 210, as shown in Figure IT51 At several times, the 
turbulence seems ready to disappear, only to spread out again. In order to confirm this surprising result, wc have 
carried out longer simulations at each of these low values of Re. Turbulence persisted over T = 4200 (in the usual 
advective time units) at Re = 220, over T = 3000 at Re = 225, and even over T = 15000 for Re = 230. 

Experiments |7| and numerical simulations in large domains (L L z = 128 x 2 x 64) |l6l | and numerical 

simulations in periodic minimal flow units (L x x L y x L z = 4x2x6) 01 have produced long-lived turbulence only for 
Re > 300. A number of studies IsL l2ll l22j have examined turbulent lifetimes as a function of initial perturbation 
amplitude, Reynolds number, and quenching rate (rate of Reynolds number decrease) in minimal flow units. In these 
studies, turbulence with a lifetime greater than T = 2000 was counted as sustained; experiments H, however, are 
carried out on timcscales several orders of magnitude longer than this. ISchmiegel and Eckhardtl ()2000|) studied the 
effect of quenching rate on turbulent lifetimes. For rates of Reynolds number decrease comparable to ours, they found 
that turbulence could in some cases subsist to Re = 280 or 290 for times on the order of T = 1000 to 10000; for 
quenching rates ten times faster than ours, turbulence was occasionally sustained to Re = 240. 

If we compare our results to the previous simulations, then the conclusion would be that turbulence is favored by 
a short streamwise direction L x = 10 and a long spanwise direction L z = 120. When either of these two conditions 
are lifted, the turbulence disappears. We note, however, that our simulations do not systematically vary the initial 
conditions and thus do not determine the probability of long-lived turbulence at these low Reynolds numbers near 
220. 

We also note that Toh and colleagues [23L |2^ | have recently reported results from simulations of Couette flow in 
domains with long spanwise extent compared with the MFU geometry. These simulations are for higher values of Re 
than those considered here. 

We observe an approximately periodic oscillation in time, shown on the right of Figure IT31 The oscillation period 
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FIG. 15: Evolution for 9 = 0° with L x i = L m = 10 and L z i = L z = 120. Re is lowered by steps of ARe = 5; not all intermediate 
values of Re are shown. The flow continues to have turbulent regions far below Re = 300. Left: evolution over the entire 
domain, showing the formation and disappearance of turbulent domains. Right: evolution of w(x — 0, y = 0, z = 60, t), showing 
irregular periodic cycles. 

of about 200 time units has the same order of magnitude as the minimum turbulent cycle |10( , but further analysis 
of our results is required before we can identify the streak and streamwise- vortex cycle in our flow. 

VI. SUMMARY 

We have used an extension of the minimal-flow-unit methodology to study large-scale turbulent-laminar patterns 
formed in plane Couette flow. Turbulent-laminar patterns are obtained as solutions to the Navier-Stokcs equations 
in domains with a single long direction. The other dimensions are just large enough to resolve the inter-plate distance 
and to contain an integer number of longitudinal vortex pairs or streaks. We have presented various visualizations of 
the computed turbulent-laminar patterns as well as space-time plots illustrating the formation and dynamics of these 
patterns. The time-averaged modulus of the spatial Fourier spectrum is shown to provide a quantitative diagnosis of 
the patterns. Periodic, localized, and intermittent states occur in our simulations where similar states are observed 
experimentally. 

We have explored the patterns' dependence on Reynolds number, domain length and tilt angle. The patterned 
states do not appear to depend sensitively on how the turbulence is initialized nor on the route taken to a particular 
point in parameter space. It is, however, possible that some parameter combinations may support different numbers 
of turbulent bands (although we have not yet observed this) . All states are bistable with respect to laminar Couette 
flow and if parameters are changed too abruptly, then reversion to laminar Couette flow occurs. 

It appears that large-scale patterns are inevitable intermediate states on the route from turbulent to laminar flow 
in large aspect-ratio Couette flow. A key open question is what mechanism causes laminar-turbulent patterns. These 
patterns are not only interesting in and of themselves, but may provide clues to the transition to turbulence in plane 
Couette flow. 
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